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The Liouville equation differs from the von Neumann equation 'only' by a characteris- 
tic superoperator. We demonstrate this for Hamiltonian dynamics, in general, and for 
the Jaynes-Cummings model, in particular. - Employing superspace (instead of Hilbert 
space), we describe time evolution of density matrices in terms of path integrals which 
are formally identical for quantum and classical mechanics. They only differ by the in- 
teraction contributing to the action. This allows to import tools developed for Feynman 
path integrals, in order to deal with superoperators instead of quantum mechanical com- 
mutators in real time evolution. Perturbation theory is derived. Besides applications in 
' classical statistical physics, the "classical path integral" and the parallel study of classi- 

C ~ 3 ' cal and quantum evolution indicate new aspects of (dynamically assisted) entanglement 

C " ^ , (generation). Our findings suggest to distinguish intra- from inter- space entanglement. 



Keywords: path integral; Liouville equation; von Neumann equation; superoperator; en- 
tanglement; Jaynes-Cummings model. 



' 1. Introduction 



The quantum-classical divide has been inte nsely studied in recent years with pro- 
found impact on various areas of research l^l^^'^i In particular, this concerns the 
foundations of quantum mechanics, new quantum technologies (quantum informa- 
tion processing, precision measurements, designer materials, etc.), recent observa- 
tions of quantum coherent processes in biology, and, last not least, unresolved issues 
surrounding "quantum gravity" E|. 

Not surprisingly, these modern topics, which touch the foundations of quantum 
mechanics in one way or another, increase the impetus to try to reconstruct and 



''I.e., the conflict between quantum mechanics necessitating an external time and diffeomorphism 
invariance in general relativity, for example, which defies its existence. Despite its great successes 
in describing the statistical aspects of experiments, quantum theory itself presents problems of 
interpretation, which are brought to the forefront in quantum cosmology. They arise from its 
indeterministic features and are clearly seen, for example, in the measurement problem. 
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2 H.-T. Elze et al. 

to better understand the emergence of quantum mechanics from simpler dynamical 
structures beneath or more profound theoretical principles. 

Indeed, there is a growing number of deterministic models of quantum mechan- 
ical objects which are based on co n.iectured fundam ental in format ion loss, coarse 
graining, or dissipation mechanisms I '^ l ^ l ^ l ^ l ^ l ^ l '*^^ ^^^, see Refs.^^^^l^for most recent 
arguments. - We recall that 't Hooft's existence theorem f^^lgjiQ-^yg thai the evolution 
of all quantum mechanical objects that are characterized by a finite dimensional 
Hilbert space can be captured by a dissipative process. This holds also for objects 
that are described by a set of mutually commuting Hermitean operators ^1^. 

However, a theory is lacking that would generally explain the emergence of 
quantum features of common objects, at the scales where they are observed. 

In order to make progress in these matters, it may be useful to further examine 
the quantum-classical divide. Presently, we look more carefully into the common as 
well as the distinctive features of classical and quantum dynamics, as described by 
the Liouville and the von Neumann equations, respectively. 

We will derive a new path integral representation of the propagator for density 
matrices in the classical theory. It is identical with the usual one a t th e kinematic 
level, employing the Feynman propagator of quantum mechanics this allows 
'external sources' in the relevant action that are coupled to terms linear or quadratic 
in the generalized coordinates. Yet the interaction part differs in a characteristic 
way0. The new formalism based on superoperators will be presented and illustrated 
here by perturbation theory applied to an anharmonic oscillator. 

Similarly, as a case study, we will re-derive the Jaynes-Cummings model 23 - 
the well-known benchmark model of quantum optics and cavity QED - based on 
classical dynamics described by a Liouville equation. Thus, when applied to the 
two-level dynamics of Rydberg atoms coupled to one mode of the photon field 
in a suitably tuned cavity, we find surprisingly that it is "almost classical" , with 
quantum and classical dynamics differing by a characteristic superoperator. 

We conclude by pointing out some interesting problems, concerning the prepa- 
ration of entangled states, in particular. Here, the quantum-classical divide shows 
new aspects, which may help to further unravel the hidden dynamics that must be 
involved when it is crossed - be it in the "classical limit" or following axiomatic 
"quantization rules" . 



2. Hamiltonian dynamics and the Liouville superoperator 

To begin with, we will consider an object with a single continuous degree of freedom. 
We will treat an atom interacting with the electromagnetic field in a later chapter, 
while a relativistic field theory has been studied elsewhere 1^^. 

'^Another path integral for classical mechanics exists, which implements Hamilton's equations as 
constraints, see Refs.Eland references there. In this approach, an analogy with quantum mechanics 
consists in the path integral as such, yet its integrand bears no resemblance. 
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A path integral for classical dynamics, entanglement, and Jaynes-Cummings model at ... 3 

Let us assume that there are only conservative forces and that Hamilton's equa- 
tions are determined by the generic Hamiltonian function: 

H{x,p) :=lp2^y(a:) , (1) 

defined in terms of generalized coordinate x and momentum p (a mass parameter 
will be inserted in due time, but is omitted here for simplicity), and where V{x) 
denotes the potential. - An ensemble of such objects, for example, following trajec- 
tories with different initial conditions, is described by a distribution function p in 
phase space, i.e., by the probability p(x,p; t)dx(ip to find a member of the ensemble 
in an infinitesimal volume at point {x,p). This distribution evolves according to the 
Liouville equation: 

with V'{x) := dV{x)/dx. - We recall that the relative minus sign in the Poisson 
bracket, or between terms here, reflects the symplectic phase space symmetry. This 
will translate into the familiar commutator structure in Eq. ([5]). 

A Fourier transformation, p{x,p; t) ^ / dy e~^Py p{x, y; t), replaces the Liouville 
equation by: 

idtp= {-dyd^+yV'{x)}p , (3) 

without changing the symbol for the distribution function, whenever changing vari- 
ables. Thus, momentum is eliminated in favour of doubling the number of coordi- 
nates. Finally, with the transformation: 

Q:^x + y/2 , q:^x-y/2 , (4) 

we obtain the Liouville equation in the form: 

idtp ^ {HQ-Hg+£{Q,q)}p , (5) 

H^-.^-^d^ + Vix) , for x = Q,q , (6) 

£{Q, q) {Q - q)V'{^^) ~ ViQ) + V{q) - -£{q, Q) . (7) 

We remark that the presented reformulation of classical dynamics is rather inde- 
pendent of the number of degrees of freedom. It applies to matrix valued as well 
as to Grassmann valued variables, representing the "pseudoclassical" fermion fields 
introduced by Casalbuoni an d by Berezin and Marinov. Field theories require a 
classical functional formalism ISCH. 

Furthermore, the Eq. ([S]) appears as the von Neumann equation for a density 
operator /5(i), considering p(Q,q;t) as its matrix elements. We automatically re- 
cover the Hamiltonian operator H related to the Hamiltonian function, Eq. ([T]), as 
in quantum theory. However, an essential difference consists in the interaction £ 
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4 H.-T. Elze et al. 

between bra- and ket- states. The Hilbert space and its dual here are coupled by a 
superoperator^ 

Since the interaction £ is antisymmetric under Q <r^ q, the complete (Liouville) 
operator on the right-hand side of Eq. ([5]) has a symmetric spectrum with respect to 
zero and, in general, will not be bounded below. Therefore, with this coupling of the 
Hilbert space and its dual by the superoperator, corresponding to the absence of a 
stable ground state, our reformulation of Hamiltonian dynamics does not qualify as 
a quantum theory. Related observations were discussed, for example, in Refs.' ^''^ ' ^ - ^! 

However, the following fact has been discussed in Refs.^^^l 

£ = potential V{x) is constant, linear, or harmonic , (8) 

with an eye on the possibility of having quantum phenomena emerge due to discrete 
spacetime structure. Analogously, the vanishing of £ in a field theory is equivalent 
with having massive or massless free fields, with or without external sources, and 
with or without bilinear couplings. Generally, in these cases, anharmonic forces or 
interactions are absent. 

In the following main parts of this work, we will study in more detail the classical 
Hamiltonian dynamics described by Eq. ([S|), or by its appropriate generalizations, 
and pay particular attention to the presence of the superoperator £, when compar- 
ing with the von Neumann equation. 

Concluding this introductory section, we recall relevant aspects of the interpre- 
tation of the density operator p{t), which we invoked here. 

2.1. Expectations, operators and the Born rule 

We begin with the normalization of the classical probability distribution: 

1 = y ^ p{x,p; t)= I dQdq S{Q - q)p{Q, q- 1) -: Tr p(t) , (9) 

incorporating the transformations of Section 2. Consider a complete set of or- 
thonormal eigenfunctions of the operator of Eq. ([51), defined by gj(x;t) := 
exp{—iEjt)gj{x) and H^gj{x) = Ejgj{x), respectively, with a discrete spectrum, 
for simplicity. Then, we may expand p: 

p{Q,q;t) = J2P3f^(*^9AQ;t)gl{q;t) . (10) 
Employing this, the normalization condition (|9]) can be stated as: 

1 = E P^tK'^'"''''''^' J dO 9j{Q)9k{Q) = E Pnit) ■ (11) 

'^This superoperator is of a very specific form, which leads to the antisymmetry in Eq. jTjl. It differs 
from a Lindblad superoperator, often obtaine d as a symmetric double commutator structure, in 
the case of open quantum mechanical systems 1121. 
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A path integral for classical dynamics, entanglement, and Jaynes-Cummings model at ... 5 

Since the classical phase space distribution is real, the expansion coefficients form 
a Hermitean matrix, pij = which we also denote by p. 
The classical expectation values are calculated as follows: 

"^xpix^p-.t) = J dQdqS{Q-q)^^^p{Q,q-t) , (12) 
= :Tr {Xp{t)) , (13) 



(p) ■■= J -^pp{x,p-t) = j AQAq 5{Q - q){-i) '^^ V (Q,g;^) , (14) 

=:Tr(Pp(t)) , (15) 
introducing the operators X and P, with matrix elements X{q^ Q) — 5{Q — q){Q + 
q)/2 and P{q,Q) = —i{5{Q — q) dq — dq S{Q — q)) (derivatives act left or right, 
as indicated). Eliminating one of the two integrations in the above equations with 
the help of the 5-functions and suitable partial integrations, these operators are 
recognized as the coordinate and momentum operators of quantum theory. 
Similarly, we find: 

^xpp{x,p;t)^^Tr{{XP + PX)p{t)) , (16) 

which constitutes an example of the symmetric Weyl ordering, when replacing clas- 
sical phase space quantities by quantum operators. - However, we remark that 
Hilbert space operators appear here by rewriting classical statistical formulae and 
not by following a quantization rule. 

The Eqs. ([T^ - (|16l) are in accordance with the interpretation of p{Q,q;t) 
as matrix elements of a density operator p(t). - However, there is an important 
caveat: The eigenvalues of normalized quantum mechanical density operators are 
usually constrained to lie between zero and one, corresponding to the interpretation 
as standard probabilities. This is not necessarily the case with the operator p ob- 
tained from a classical probability distribution. Similarly, the Wigner distribution 
- obtained from the matrix elements of a quantum mechanical density operator by 
applying the transformations leading from p(x,p) to p{Q,q) in reverse - generally, 
is not positive semi-definite on phase space, even though its marginal distributions 
are. Therefore, it does not qualify as a classical probability density. 

As we have indicated before, there is clearly a dynamical feature missing, which 
governs the crossing of the quantum-classical divide, if not done 'by hand', as in 
any of the usual "quantization prescriptions". Last not least, this must establish the 
Born rule by eliminating negative probabilities or by leading to their satisfactory 
interpretation. 

3. From Hilbert space to superspace 

In this section, we reformulate the notions relevant for the dynamics of density 
operators, at which we arrived in the previous section, in a more convenient way. 
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6 H.-T. Elze et al. 

introducing the concept of superspace El, also called Liouville space - see Ref.'^ for 
a concise presentation and numerous applications. 

Considering a physical object characterized by the Hamihonian H , as in quan- 
tum theory, we introduce a complete set of basis states, {|j)} (j — 1,...,A^), 
assuming that the relevant Hilbert space is iV-dimensional. Then, taking matrix 
elements of the von Neumann equation, for example, we have: 

idtpjk = [[Hp)jk ~ [pH)jk] , j, fc = 1, 2, . . . , iV , (17) 

with a density matrix p of iV^ elements. Which may be written as: 

idtPjk —''^^Cjk,lmPlm , (18) 
l.m 

in terms of the Liouville superoperator C defined by: 

^jk,lm ■= HjiSkm — Hl^5jl . (19) 

This definition suggests to introduce a space where the density operator is a vector, 
which is the Liouville space (or superspace). The dynamics of the density operator 
can then be more conveniently described, completely in parallel for classical and 
quantum mechanics, as we shall see. 

Given the Hilbert space, as above, the density operator can be expanded as: 

P = Y.P^^\^)^^\ ■ (20) 

We may think of the family of N'^ operators \j){k\, with j, fc = 1, • • • , TV, as a 
complete set of matrices, or vectors, and express the density operator as: 

|P>=E'^^'=I^'^)) ' (21) 

where the "ket" \jk)) denotes the Liouville space vector representing the Hilbert 
space operator |j)(fc|. Similarly, we introduce a "bra" vector {{jk\ as the Hermitean 
conjugate of \jk)). 

In Liouville space, any operator A is represented by a vector and denoted by 
\A)). It can be expanded as: 

\A))=Y,A,,\jk)) , (22) 

3,k 

where Ajk are the usual matrix elements (j|^|fc). - Furthermore, we can define a 
"bra" vector {{B\, representing , and the scalar product of two operators: 

{{B\A)) := Tr(BU) . (23) 

'^The notion of superspace here, at first sigfit, lias little in common and sfiould not be confused 
with the corresponding term relating to supersymmetry. 
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A path integral for classical dynamics, entanglement, and Jaynes-Cummings model at ... 7 

Then, the following orthonormahty condition holds: 

{{]k\-mn)) = Tr(|A:)(j|m)(^l) = 4.<5jm , (24) 
which is analogous to {j\k) — Sjk- - Finally, consider the scalar product: 

{{jk\A)) = Tr(|fc)o-|i) = Y^muim = oiiifc) ^ a,, . (25) 

Upon substitution in Eq. (f^ . this yields: 

\A)) ^J2\^k))mA)) . (26) 

This is consistent with the following completeness relation in Liouville space: 

J2\jk)){{jk\^l . (27) 

Following these considerations, it can be verified that Liouville space is a linear 
space, in which the density operator p is a vector. In this space, a linear operator 
can be defined by: 

J':= \jk)){{jk\T\mn)){{mn\= ^ Tjk.rnn\jk)){{mn\ , (28) 

j,k,m.n j.k.m,n 

i.e., in terms of appropriate matrix elements. 

The importance of Liouville space for classical and quantum dynamics is that 
the Liouville and von Neumann equations, both, can be written in the form: 

idtp^tp , (29) 

with an appropriate superoperator C, cf. Eqs. ©-(O and ([T7l) -(fT9 |) . respectively. 
Thus, there is a formal analogy (even isomorphism) between the structure of these 
equations and the Schrodinger equation, idt'^ = H'i' . 

These observations suggest that techniques or formal results concerning the 
solution of the Schrodinger equation can be transferred to the case of the Liouville 
or von Neumann equations with the help of Liouville space notions. This concerns 
perturbation theory (and nonperturbative methods) as much as a path integral 
formulation, which we shall discuss in turn. 

First of all, we introduce the Liouville space evolution operator U satisfying: 

tdtU{t,to) = mU{t,to) , (30) 
with the initial condition U{t,to) — 1. This implies: 

m^^it,to)p{to) , (31) 

as the solution of the density operator equation of motion. For a time independent 
Liouville superoperator this yields: 

U{t,to)^exp{^iC{t~to)) . (32) 
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8 H.-T. Elze et al. 

Thus, time evolution of the density matrix is implemented by a superoperator in 
Liouville space, while in Hilbert space the evolution is described by: 

pit) = U{t,to)pito)UHt,to) , (33) 

with U{t,to) := exp{—iHt). The Eqs. (I5T]) and represent the evolution of the 
same object, although in different spaces. - For a time dependent Hamiltonian, we 
have instead: 

U{t,to)^Texp{~i [ drCir)) (34) 



- 1 + / dT„ / dT„_i . . . / dn £(t„)£(t„_i) . . .£(ti) ,(35) 

n=l "^'o •'to Jto 

in terms of the time-ordered exponential. 

For later purposes, we also introduce the "interaction picture" in Liouville space. 
Considering a Liouville operator which consists of two parts: 

£ = £o{t) + C\t) , (36) 

we obtain the evolution operator in the following form: 

Uit,to)=Uoit,to)Uiit,to) , (37) 

with: 

4(i,^o) = Texp(-i /" dr£o(T)) , (38) 



and: 

Ui{t,to)^Texp{-i f drC'jiT)) , (39) 

Jto 

with C'j{t) := UQ{T,to)£'{T)Uo{T,tQ). For an operator Uq that can be treated ex- 
actly, study of time evolution essentially concerns the operator Uj - for this, Eq. (|39|) 
presents the starting point of perturbation theory (expanding the exponential). 

4. The quantum/classical path integral for the propagator of 
density matrices 

The technical ingredients needed for the Feynman path integral approach, for the 
derivation of quantum mechanical propagators in particular, are very well known. 
We import these here, especially from Ref.'^, in order to derive a path integral for 
the propagator of density matrices based on the Liouville space formulation of the 
preceding Sections. 

Our derivation relies on the close formal similarity between the classical Liouville 
equation and the von Neumann equation on one hand side and the Schrodinger 
equation on the other, in an appropriate representation, as we have discussed 0. 

'^In this chapter, we reinstate h explicitly. 
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A path integral for classical dynamics, entanglement, and Jaynes-Cummings model at ... 9 

4.1. Essentials of the Feynman path integral 

We recall that the (forward propagating) operator Green's function G, 

Git,to) = eit-to)ex_p{-iH{t^to)/h) , (40) 

allows one to write the solution of the time dependent Schrodinger equation as: 
\'tp{t)) = G{t, to)\^{to)). - Correspondingly, for t > to, we have the coordinate space 
matrix elements: 

G{x,t;y,to) ^ {x\e-^^^'-'°^/''\y) , (41) 

from which one derives the path integral representation of these amplitudes, for 
a generic Hamiltonian H = fp' /2m + V(x), describing a pa rtic le of mass m in an 
external potential V , cf. Eq.(IT]), through the following steps 

• Cut the time interval from t^ to t into a large number N of equal pieces. 

• Write the exponential of the Hamiltonian operator x time as a product of 
identical factors, each factor representing the propagator for a small time 
interval (xl/N . 

• Separate the kinetic and potential terms contributing to H in each one of 
these factors with the crucial help of the Trotter product formula. 

• Alternatingly, insert complete sets of momentum and coordinate eigen- 
states, such as J dx \x){x\ = 1 (and correspondingly for momentum eigen- 
states) between the factors of exponentials involving either momentum or 
coordinate operators and evaluate the resulting Gaussian integrals over 
momentum variables. 

• Realize that the obtained phases in the product of exponentials can be 
summed up to represent a discretized version of the classical action per- 
taining to the Hamiltonian function (corresponding to H). 

Taking the limit — )■ cx) in the end, one obtains the following Feynman path 
integral representation of the amplitudes in question: 

G{x,t-yM) = (42) 

. . N-l 

lim i-^r'^ / dxi...dx^.i exp(f {^C-'^^^f-yi-^))) 

Vx exp Q / dr - y(x))) (43) 

= J Vx exp{^S[x,x]) , (44) 

where e := {t — toj/N, x := dx/Ar and where S denotes the relevant classical action, 
which is to be evaluated for each one of the paths contributing to the integral, with 
the boundary conditions xit) = x and x(to) = U- 
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10 H.-T. Elze et al. 

4.2. The Liouville space propagator as a path integral 

We are now ready to appreciate the economy of the Liouville space representation 
introduced in Section 3. In particular, the formal solution of the classical Liouville 
equation as well as of the quantum mechanical von Neumann equation, both, can 
be written in the form: 

Wm = e-'^('~*")/"|p(io)» , (45) 

where H is the appropriate super-Hamiltonian. Generally, we have: 

{{Q,q\n\Q',q'))^S{Q~Q')6{q~q'){H{Q)~H{q)+£{Q,q)) , (46) 

where H denotes the appropriate Hamilton operator in coordinate representation, 
as indicated, which alone is relevant for the von Neumann equation, while £ repre- 
sents the additional superoperator for classical dynamics, cf. Section 2. 

In order to solve the problem of time evolution in the present case, we need to 
know the (super)matrix elements entering the propagation equation: 

{{Q,q\pm - / dQ'dq' ((Q,g|e-^«(*-*«)/'"'|Q',q')>((Q','z'|p(io))) , (47) 

which appears formally analogous to evolution of a state vector according the 
Schrodinger equation. Thus, not surprisingly, we go through the steps indicated 
in the preceding Section 4.1, in order to construct the path integral representation 
of the propagator here. 

However, in this derivation, we have to pay attention to the crucial role of the 
Trotter product formula. It turns out that it can be generalized for our purposes, 
where superoperators present the new feature, in a straightf orw ard way; the relevant 
definitions and details of the proof will be given elsewhere 

Rewriting the Eq. dJ?]) as: 

P(Q, q;t)^ J dQ'dq' GiQ, q- t\Q' , q'; to)p{Q\ q'; to) , (48) 

our interest is to know the superpropagator G- Next, we will follow the recipe to 
arrive at a path integral representation, as summarized above, in Section 4.1. 

In particular, here we make use of suitably inserted complete sets of superspace 
vectors, such as: 

I dQdq\Q,q)){{Q,q\ = l , (49) 

and, correspondingly, for momentum space, cf. Eq. ([27]) ■ Using the plane wave re- 
lation between coordinate and momentum eigenfunctions, we also employ: 

mp\ = ^/ dQd<7 exp(- ^(PQ-m))((Q,'?I • (50) 

Furthermore, the orthogonality relation {{Q , q\Q' , q')) = 5{Q — Q')5{q — q'), cf. 
Eq. dMI, implies: 

{{P,p\Q,q)) = ^exp{-^iPQ-pq)) . (51) 
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A path integral for classical dynamics, entanglement, and Jaynes-Cummings model at ... 11 

Then, with all following steps of the derivation in parallel with the u sual ones in 
quantum mechanics, it is straightforward to obtain in the present case^^^l 

g(Q/,g/;t|Q„gi;to) = J VQVq exp {^S[Q,Q;q,q]) , (52) 

with the boundary conditions Q(ti) — Qi, q{ti) = g^, Q{tf) = Qf, and q{tf) = qf, 
and where the superaction 5, corresponding to the super-Hamiltonian %' above, is 
defined as follows: 

S = f dr {nQ,q)-V{Q,q)) (53) 

J to 
J to 

We recall that £ = corresponds to evolution according to the von Neumann equa- 
tion, whereas £ 7^ represents classical dynamics in accordance with the Liouville 
equation, cf. Eqs. ©-Q. 

However simple this result may seem, the Eqs. (j52l) ~ (j54p present a new approach 
to describe time evolution of the full density matrix, with the particular feature that 
classical and quantum mechanical motion are formally treated in parallel, differing 
only in the action entering the phase in the path integral]^. 

We emphasize that our derivations are not confined to one-dimensional or single- 
particle physics, but can be extended as well all the way to relativistic field theories. 
Various applications come to mind here, some of which will be discussed in the 
following and in the concluding section. 



4.3. Perturbation theory and superpropagator Dyson equation 

Considering the splitting of the superaction as in Eq. (|53p . the perturbation the- 
ory naturally departs from organizing contributions to the full superpropagator, 
Eq. ([52]) . according to powers of the "perturbation" V. Sometimes it may be ad- 
vantageous to include parts of the perturbation into the "free" part T. This must 
be familiar from quantum mechanics, which presents a special case of our general 
considerations here. 

To begin with, if V{Q,q) = V{Q) — V{q), corresponding to the superoperator 
related to the von Neumann equation, then the path integrals in Eq. (|52p factorize 
and we recover quantum mechanics. 

In the absence of an external potential or other interactions (V = 0), the zeroth 



'We remark that a complementary approach, based on the effective action generating equ al-time 
correlation functions for nonequilibrium statistical systems, has been presented in Ref.l^^l^ which 
results in evolution equations for a truncated set of correlation functions. 
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12 H.-T. Elze et al. 

order or free superpropagator Qq is obtained as: 



Go{Qf,qf;t\Q^,q^;to) = J VQVq ^^"Pi^j^ dTT{Q,q)) (55) 
= Go{Qf,t;Qi,to)Go{qf,t;q.i,to) , (56) 



in terms of the well known free quan tum mechanical propagator Go, cf. Eqs. (|42p - 
((44|) , which is explicitly given by 

Go{x,t;y,to)^Go{x,y;T:=t-to) = {^^y^\^p[^{x-yf) , (57) 

for a free nonrclativistic particle of mass m. Remarkably, this zeroth order result is 
identical for classical and quantum mechanical propagation. 

Despite the fact that the free propagator for the Schrodinger equation incor- 
porates such phenomena as the quantum mechanical spreading of a wave packet, 
we learn here that it also describes the propagation of a free classical particle. It 
is straightforward to verify - following the transformations between Eqs. ([T]) and 
Eq.([7|) - that a massive particle, initialized as p{x,p, to) :— 2ttS{x — xo)6{p — po) is 
propagated to p{x,p,t) :— 2'!tS{x — xq — Tp/m)5{p — pq), as expected. 

For the perturbative expansion, we employ the standard formula: 



exp ( - ^ /* dr V(Q, q)r) =f]-A^ V(Q, q)r 

" Jo „_n ^' " "'0 



(58) 



= E(t) / driV(Q,g),, ... / dT„V(Q,q),„ , (59) 

n=0 ''^ 

with V(g,g),, :-V(g(Tfc),g(Tfc)). 

In order to analyze such terms at a given order, we make use of the impor- 
tant semigro up property of the (free) propagator and obtain to first order in the 
perturbation ^\ 

Q{Q, q; t\Q', q'; to) = Go{Q, q; t\Q', q; to) 



to 



dr / dxdy GoiQ, q; t\x, y; t)V{x, y)go{x, y; t\Q' , q'; to) + 0{V^) , (60) 



to be illustrated explicitly by the result for an anharmonic potential shortly. 

We remark that on the right-hand side of Eq. (|60|) the superpotential is preceded 
(and followed) by a zeroth order propagator. This observation, which similarly holds 
at every order of this expansion, leads to a recursion relation of the fc-th order prop- 
agator in terms of the (fc— l)-th order one. This allows us to resum the perturbation 
series in the form of a Dyson integral equation for the full superpropagator: 

5iQ, q; t\Q', q'; to) - GoiQ, g; t\Q', <?'; to) 

~'^It J ^^'^y ^o{Q,q;t\x,y;T)V{x,y)g{x,y;T\Q' ,q';to) ■ (61) 
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The whole procedure follows the usual one in quantum mechanics, yet includes the 
case of classical mechanics, and possibly others, for a suitably chosen superaction. 

4.4. Illustration: the case of an anharmonic potential 

In order to make our general derivations more concrete and to extract some inter- 
esting general aspects, it may be useful to consider the example of a massive particle 
in an anharmonic potential, V{x) :— Xx*, where A is the coupling constant. - We 
recall that for constant, linear, or harmonic coupling terms there is no difference 
between classical and quantum dynamics, cf. in the representation that we have 
developed in this article. 

The calculations evaluating the superpropagator to first order, here with: 

Vix,y)^{x-y)V'{^)^^{x^-y^ + 2{x'y-xy')) , (62) 

for classical dynamics {V{x,y) = V{x) — V{y) for quantum mechanics) consist in 
straightforward (if tedious) multiple Gaussian integrals, according to Eqs. ([55)) -((57 )) 
and Eq. IpO]) . The final result is: 

g{Q,q;t\Q',q';to) = QoiQ , q;t\Q' , q' ■,to) 

\ - jX[CiTqM{Q, q- Q\ q'- T) + C^T^^i^iQ, q; Q' , q'; T)] ) + 0{X^) , (63) 

where T := t — to, the coefRcients for classical dynamics, Ci := 1/2, C2 := 1/2 
(Ci :~ 1, C2 for quantum mechanics), and with the function: 

rQM(Q- 9; Q', 9'; T) := ^ [i — (3Q2 + 4QQ' + 3Q'') 



T 



+Q^ + Q'Q' + Q^Q'' + QQ" + Q'--{Q,Q')^{q,q') , (64) 

J 5 L J 
where the term is repeated, as indicated, with an exchange of variables and complex 
conjugation. Similarly: 

T / ihT , 
5 V m 



rcL(Q' 1; Q', q'; ^(3^9 + 2Qg' + 2Q'q + 3Q'q' 



+ i [Q3(4g + q') + Q^Q\3q + 2q') + QQ'\2q + 3q') + Q'\q + Aq')] 

-\[{Q.Q')^{q.q')\ ) . (65) 

This result shows several interesting features. First of all, the perturbative ex- 
pansion turns out to be a short-time expansion, with the overall scale of the first 
order correction set by AT. Furthermore, different contributing terms differ by a 
scale set by T/to, i.e., by T x Compton wavelength of the particle. Numerical 
studies visualizing the outcome here are presently underway 

However, most interesting seem general similarities and differences between clas- 
sical ( "CL" ) and quantum mechanical ( "QM" ) result in Eq. . The CL result has 
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the same zeroth order term as QM; at first order, CL has one term in common with 
QM which, however, is reduced by an overall factor 1/2. This obviously stems from 
the varied expressions for V between CL and QM, cf. Eq. (|62|) . For the same reason, 
CL has additional terms, collected in Fq]^, Eq. (|65|) . which are absent in QM. 

4.5. Intra- and inter-space entanglement 

There is a qualitative difference between CL and QM, contained in Fqj^ and based 
on the different superoperators that enter the full path integral, Eqs. (f52|) - (|54|) . 
Equivalently, since the QM evolution is generated by a commutator of the Hamil- 
tonian with the density operator p, it superposes and, for multi-partite systems, 
generally, entangles underlying bra- and ket-states separately, cx Hijpjk — pijHj).. 
For a bi-partite system, it is revealing to write such terms more clearly as: 

[Hint,p] = Hipi (g) H2P2 - PiHi (8) P2H2 , (66) 

for an interaction (x Hi ^ H2, with the factors acting on subsystems "1" and "2", 
respectively, and where p — pi®p2, for a separable initial state. This ha s been ca lled 
dynamically assisted entanglement generation, see, for example, Refs. I 23 | 24 | 25 ] 

It may come as a surprise that the CL evolution does this just as well, due to 
the contribution of Fqjyj for the first two terms on the right-hand side of Eq. (|62|) 
or, generally, due to the superoperator £ of our earlier considerations. For polyno- 
mial interactions, for example, this superoperator always contains a contribution 
proportional to the usual QM terms. 

However, the CL evolution produces additional correlations in p, due to the 
generator cx Cij-kiPki, which possibly entangles bra- and ket-states. - In comparison 
with Eq. (|66p . for example, such terms can have the unfamiliar structure: 

H[pi (E) P2H2 - PiH[ ® H!iP2 , (67) 

which differs decidedly from a commutator. - This leads us to distinguish intra- 
(i.e., within given tensor product Hilbert space of subsystems "1" and "2") and 
inter-space entanglement (i.e., between said Hilbert space and its dual). 

For example, consider the anharmonic potential V{xi — X2) '■= \{xi — X2)'* 
for a bi-partite system consisting of particles "1" and "2" . Following and suitably 
generalizing our derivation in Section 2, this leads to the interaction: 

V(Qi,Q2;gi,<Z2)-^A(Qi-gi-(Q2-g2))(0i + gi-(Q2 + g2))' , (68) 

in terms of variables introduced previously, taking into account both subsystems; 
similarly as before, the Q and q variables refer to bra- and ket-states, respectively. 
Besides the separable terms, oc {Qa — qa){Qa + QaY , o = 1, 2, there are the terms 
which mix (and entangle) variables of both subsystems, as usual in QM. However, 
there are clearly additional terms that refer to Hilbert space and its dual simulta- 
neously (and entangle corresponding states), for example, esc QaQbQb^ b ^ a. 
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In retrospect, somehow, such difference between CL and QM evohition had to 
be expected: instead with superstates \Q,q)), we could have worked with super- 
states |a;,p)), relating to coordinates and momenta of the classical theory. There, 
coordinates and momenta end up tightly correlated, due to Hamilton's equations, 
and produce inter-space entanglement in an interacting bi-partite system. 

Thus, we find that the confrontation of CL with QM, as in our side-by-side 
study, is quite revealing. In particular, we speculate that this opens new views 
on generating entanglement in multipartite systems, perhaps, by evolving through 
quasiclassical stages or by making use of decohered intermediary states [3. 

Concerning the quantum-classical divide, the present analysis shows that there 
is a deep formal similarity between CL and QM. However, this also demonstrates 
that what has been discussed in various ways as CL limit of QM - and which is 
similarly relevant for "emergent QM" - deserves more study. 

While our work has been concerned mainly with the evolution of CL or QM 
objects, we recall that V.I. Man'ko and collaborators have pointed out that classical 
states may differ widely from what could be obtained as the "?i — >■ 0" limit of 
quantum mechanical ones. They show that all states can be classified by their 
'tomograms' as either CL or QM, CL and QM, and neither CL nor CM 

The classical limit might be a "For^llftacticalPurposcs" limit, gradually ap- 
proached through decoherence or "h — > 0" . However, in order to bridge (if at all) 
the qualitative difference between intra- and inter-space entanglement that we find, 
and explain the "Man'ko classes of states" , some unknown dynamics beneath still 
awaits to be uncovered R. 



5. The almost classical Jaynes-Cummings model 

In this section, we apply our operator approach for the Liouville equation to a field 
theory, namely to a Rydberg atom interacting with the electromagnetic field. Fol- 
lowing the approximations that lead to the quantum mechanical Jaynes-Cummings 
model we will show that the dynamics of this celebrated model is almost of 
classical character. As we shall see, if it were not for the anharmonic Coulomb 
interaction between electron and atom, the dynamics would be entirely classical. 



^Previous considerations of the semiclassical regime, such as in Refs. | 23|24| ^ ^g[.g motivated as 
suitable approximations of the quantum mechanical evolution, in particular, for studies of the dif- 
ferent decoherence properties between classically regular and chaotic systems. Our results seem to 
show that crossing the quantum-classical divide may offer an additional resource for entanglement 
generation and related "truly quantum" phenome na. Th is might be related to the "underlying 
reality" of (CL and QM) physics, assumed in Refs.l2Sl2ZI, consisting in statistical correlations. 
''A simple attractor model, motivated by ass umptions about effects of fundamental spacetime 
discreteness has been discussed in Ref. Ls^. 
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5.1. The classical model 



We consider an electron (mass to) interacting electromagnetically with a positive 
charge (atom) fixed at the origin and with the radiation field. Thus, we depart from 
the classical Lagrangian: 

L:=^i' + Jd'r{l-iE'^B')^p^ + J-A} , (69) 

where the electric and magnetic fields, respectively, are given by: 

E ^ -A-V(l) , B:^V xA , (70) 

as usual, in terms of vector and scalar potential, A and 0, respectively. The charge 
and current densities, p and J, respectively, are given by: 

p{r) :=^-ed^{r-x)+6^{r) , J{r) ^ -exS^{r - x) . (71) 

Next, we introduce Fourier modes of the fields, with A{k) = A*{—k) and (j){k) = 
(f)* {—k), since the fields are real. We choose the Coulomb gauge by imposing A|| (k) = 
0, which implies V • A = 0. Correspondingly rewriting the Lagrangian, we determine 
the canonical momenta, in order to obtain the Hamiltonian of the classical model: 

H = ^(p + e J d^k {^(fc)e''=-^ + A*(fc)e^*'=-^})^ - 

+_L J {n*(fc) • n(fc) + k^A*{k) ■ A{k)} (72) 

= H{x,p;A,n*;A*,n) , (73) 

where we indicate the canonically conjugated pairs of variables of the Hamiltonian; 
the momentum integrations have to take into account that not all Fourier modes 
are independent, for real fields. 

We are now in the position, cf. Section 2, to describe this model in phase space. 
We proceed in four steps: 

• First, we introduce the probability density (over phase space) 
p{x,p; A,Il*; A* ,IV), which will be interpreted, as before, as matrix el- 
ement of a Hermitean density operator p. We assume that the atom- 
electromagnetic-field system is confined to a cavity of finite volume V, 
thus replacing integrals by discrete mode sums, / d^k g{k) — > V^^ gt- 

• Second, we obtain the Liouville equation, —dtp — {H, p} ~ . . . , evaluating 
the relevant Poisson bracket. 

• Third, we replace momenta by coordinates via Fourier transformation, p — >■ 
y, n* ^ Bk, Uk ^ B*. 

• Fourth, we perform the "Wigner rotations", Q := x + y/2, q := x — y/2, 
Qk := Ak + Bk/2, and qk := Ak - Bk/2. 

Details and the following derivations will be reported elsewhere 
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Introducing the following notation: 

vix) ■■= , for X = Q, g ; q) ^^'^^ - viQ) + ^ii) . (74) 

where we suppress constant normahzation factors etc., as before, we obtain the 
remarkable result that the classical evolution equation is: 

idtp = { "von Neumann" + £ + V + Yj]p (75) 

-{-i^cf + nQ)-(-i5/ + n<z)) 

k 



e 

k 



_^_V[e*'?Qfe-5Q+e*'gfc.9,] 



k,k' 

+ {£{Q,q)+T + j:}p , (76) 

with p = p{Q, QkiQ*k'^ 9; qk,qk\t), where k runs over all modes, ujk :— \k\, and where 
r and S denote rather complicated terms that involve all phase space variables; they 
are given explicitly in Ref. While the last line of Eq. ([76]) presents additional 
terms, in particular the superoperator the previous terms represent exactly the 
terms of the quantum mechanical von Neumann equation for the atom- field system 
under consideration; besides further interaction terms, due to minimal coupling, we 
find the contribution of the free electromagnetic field in the second and those of the 
electron interacting with the Coulomb potential of the Rydberg atom in the first 
line, respectively. 

We anticipate that in the dipole approximation we have F, S 0. Therefore, we 
do not study further the impact of those terms here 

Instead, we recall the well known additional approximations that turn the 
von Neumann terms of Eqs. (|75p -(|76 p into those of the Jaynes-Cummings model^^^. 

• The dipole approximation, assuming that k ■ I <ti 1, where hk and / denote 
a typical photon momentum and linear size of a Rydberg electron orbit, 
respectively. 

• The restriction to one cavity photon mode with energy huj. This yields the 
approximate Hamilton operator 

H^Y. + ^{a'^a + ^) + i^d^jia - a^)\i){j\ , 

where the sums run over the Rydberg levels, with energies Hcoi, a^^^ are 
photon annihilation (creation) operators, and where the last term involves 



I 



I 



June 9, 2010 0:21 WSPC/INSTRUCTION FILE Elze-May2010 



18 H.-T. Elze et al. 

the dipole transition amplitudes dij. 

• The restriction to a two-level subspace, spanned by states \g), |e); the lower 
level energy is conveniently set to hujg = 0, while the physical realizations 
considered, usually, have w w, i.e., approximately resonant photon and 
excited electron states. 

• The rotating wave approximation, which yields the 'energy conserving' 
dipole interaction term D cx a\e){g\ — \g){e\d^ . 

The resulting Jaynes-Cummings Hamiltonian is: 

Hjc=^e\e){e\+ujid^d+^)+ideg{d\e){g\-\g){e\d^) . (77) 

Then, following the above derivation, the evolution equation becomes: 

tdtp^[HjC,P]+£p , (78) 

where we appropriately incorporated here the superoperator £. This term presents 
the only difference between the classical dynamics described by Eqs. (|75|) -(j76 |) and 
the usual quantum mechanical one. 

Thus, we find in this 'standard model' of quantum optics a detailed example for 
the similarity between CL and QM evolution laws. 

5.2. Dipole interaction and Coulomb superoperator as 
perturbations 

In order to illustrate our findings, we will briefly study the influence of the clas- 
sical superoperator on the evolution described by the Jaynes-Cummings mo del i n 
perturbation theory, while a more complete analysis will be presented in Ref. 

Following Section 4.3, cf. Eq. (|6T|) . we presently treat the dipole interaction, I) oc 
— together with the superoperator £ of Eq. fTS]) as perturbation. 

Correspondingly, we choose the Rydberg atom states and one-mode photon number 
states for the tensor product basis of the relevant Hilbert space. 

Then, the density matrix evolves according to: 

p{t) = g{t)m = go{t)m -i f dT Goit -t){V + £)go{T)m , (79) 

JQ 

to first order mV + E. While the dipole operator acts on atom and electomagnetic 
field states simultaneously, the superoperator acts only on the Rydberg states. 
Hence, the matrix elements of £ are defined by: 

£ab.cd:= J d''Q<i'qra{Q)Mq)£iQ,q)MQWd{q) , (80) 

where ipi=(n,i,m) denote standard hydrogen-like wave functions. Since these are 
eigenstates of parity. Pi = ±1, we find £ab,cd = PaPbPcPd£ab,cd, which implies 
the selection rule: £ab,cd = 0, if PaPbPcPd = —1- Furthermore, we have: Sab.cd = 

~^ba,dc- 
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Consequently, the only nonzero matrix elements are £eg,eg, £ge,ge, £ee,gg, and 
£gg,ee - The latter two vanish for the specific ground and excited states used in cavity 
QED experiments^^. 

Taking matrix elements of Eq. (|79p . we find that the superoperator only affects 
the evolution of the matrix elements Peg\nn' and pge\nn' = Ptg\n'n (^^ hermiticity) ; 
here the Fock states are labelled by photon numbers n, n' . Then, we find; 

Peg\nn'{t) = e""!--"' ([1 + zt^e^.e^lPeglnn' (0) 

+degt[y/n + Ipgg {0)p 

n+1 n' {0)^V^PeeiO)Pnn'~l{0)]) , (81) 

where we assume that the initial state factorizes, for simplicity. 

Thus, to first order in this perturbative expansion, we find that the superoper- 
ator competes with the dipole interaction, as far as the atom states are concerned; 
however, it does so without affecting the field states. Numerical estimates indicate 
that its matrix elements are not small com pared to the ones of t he dipole operator 
for cavity QED experiments (see Ref. ^^Hl for further details) 

To summarize, almost all dynamical (operator) features of the Jaynes-Cummings 
model can be derived in the classical framework, as we have shown. 

Nevertheless, there is a noticeable difference between classical and quantum 
evolution in the version of the Jaynes-Cummings model that is related to cavity 
QED experiments. This is solely due to the classical superoperator, which stems 
from the Coulomb interaction between electron and Rydberg ion. 

In distinction, had we considered a charged particle trapped by a linear or 
harmonic potential, then the superoperator would vanish identically, cf. Section 2, 
Eq. ([S]). The correspondingly modified Jaynes-Cummings model could be seen as of 
entirely classical origin, despite its quantum mechanical appearance. 



6. Conclusions 

Beginning with the Liouville equation of cla ssical statistical mechanics, we have 
introduced a (super)operator formulation which brings it as close as possible 
to the von Neumann equa tion o f quantum mechanics, provided suitable coordinates 
are chosen in superspace l ^ l ^^ i Presently, we have concentrated on the similarities 
and differences between both evolution equations . 

We have chosen the Jaynes-Cummings model l ^^ l '^^ l^ in particular, to illustrate 
both aspects and to show that this benchmark model of quantum optics a nd c avity 
QED can be interpreted to a large extent in terms of classical dynamics Fur- 
thermore, this model serves as an example that our more general considerations 
apply not only to single- or few-particle systems, but to field theories as well. 

While presently the relevant Hilbert space has been treated as tensor product 
space of the Rydberg electron single-particle and the photon Fock space, earlier also 
a functional approach combining fermion and boson fields has been discussed tl^. 

More generally, we discussed in parallel the formal solutions of the Liouville and 
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the von Neumann equations. Introducing suitable pr opagators, we derived a path 
integral representation for both cases side by side ^^U. 

The path integral for the propagator of the Liouville equation is new and may 
have interesting applications in classical physics. We derived and illustrated the 
related perturbation theory. 

We discussed how these results and the action entering this path integral, in 
particular, hint at the possibility that a form of entanglement is also generated 
by classical dynamics, which has gone unnoticed before. I t combin es the quantum 
mechanical dynamically assisted entanglement generation l ^'^ l ^^ l^ with a classical 
counterpart. We call the former intra-space entanglement, since it acts separately 
within the Hilbert spaces of bra- and ket-states. In distinction, the classical dy- 
namics additionally produces inter-space entanglement, i.e., it correlates the Hilbert 
space and its adjoint in addition to what would, otherwise, be recognized as quan- 
tum entanglement. 

If the relative strength of intra- and inter-space entanglement can be manipu- 
lated, for example, by driving a system dynamically between quantum and classical 
behaviour, this may open additional ways to influence the quantum mechanical 
entanglement, which is of central importance in research concerning quantum in- 
formation and quantum foundations alike. 

The close relation between classical and quantum mechanical dynamics that 
we uncovered may help to address in new ways problems related to the nature of 
classical or qu antuni sta tes |26 1 28|29| ^^le pathway s, if any, over the quantum- 
classical divide l ^ l ^ l ^ l ^ l ^ l, or the measurement problem 1213. 
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